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ABSTRACT 

We present new K-band spectroscopy for a sample of 48 starburst galaxies, obtained 
using UKIRT in Hawaii. This constitutes a fair sample of the most common types 
of starburst galaxies found in the nearby Universe, containing galaxies with different 
morphologies, masses, metallicities and far infrared luminosity Lm < 10 10 L@. The va- 
riety of near infrared spectral features shown by these galaxies implies different bursts 
characteristics, which suggests that we survey galaxies with different star formation 
histories or at different stages of their burst evolution. 

Using synthetic starburst models, we conclude that the best ensemble of parame- 
ters which describe starburst galaxies in the nearby universe are a constant rate of star 
formation, a Salpeter IMF with an upper mass cutoff M up = 30 M and bursts ages 
between 10 Myr and 1 Gyr. The model is fully consistent with the differences observed 
in the optical and FIR between the different types of starbursts. It suggests that HII 
galaxies have younger bursts and lower metallicities than SBNGs, while LIRGs have 
younger bursts but higher metallicities. 

Although the above solution from the synthetic starburst model is fully consistent 
with our data, it may not constitute a strong constraint on the duration of the bursts 
and the IMF. A possible alternative may be a sequence of short bursts (which may 
follow an universal IMF) over a relatively long period of time. In favour of the multiple 
burst hypothesis, we distinguish in our spectra some variations of NIR features with 
the aperture which can be interpreted as evidence that the burst regions are not 
homogeneous in space and time. We also found that the burst stellar populations are 
dominated by early-type B stars, a characteristic which seems difficult to explain with 
only one evolved burst. 

Our observations suggest that the starburst phenomenon must be a sustained or 
self-sustained phenomenon: either star formation is continuous in time or multiple 
bursts happen in sequence over a relatively long period of time. The generality of our 
observations implies that this is a characteristic of starburst galaxies in the nearby 
Universe. 
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1 INTRODUCTION 



Although a wealth of information has been accumulated on 
starburst galaxies over the last 30 years, there are still funda- 
mental questions which are left unanswered. In particular, 
we do not know if starburst galaxies form stars according 
to a universal Salpeter's initial mass function (IMF), as it 
seems to be the case in normal galaxies (Elmegreen 1999), 
or if their IMF slope is flat or has a high lower mass cutoff, 
which would allow starburst galaxies to form preferentially 
massive stars (Rieke et al. 1993). Neither do we know if the 
typical duration of a burst is short, comparable to dynami- 



cal time-scales (Mas-Hesse & Kunth 1999), or if bursts are 
spread over a longer period of time (Meurer 2000). 

In the past, various observations were carried out in or- 
der to find an answer to these two questions. Some of these 
studies have yield quite surprising results. By comparing 
different luminosity ratios in the optical and far infrared, 
Coziol (1996) found that most starburst galaxies may have 
kept forming stars at near constant rates over the last few 
billion years. The same conclusion was reached by Goldader 
et al. (1997b) for a sample of luminous infrared galaxies 
(L/H > 10 10 L ). Using synthetic starburst models, these 
authors found continuous star formation and burst ages 
varying between 10 7 and 10 9 yr. They also found that the 
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best fitted IMF is Salpeter like, but have an upper mass 
cutoff M„ p = 30 Mq . This kind of IMF in starburst galaxies 
was previously reported by Doyon, Puxley & Joseph (1992), 
which also suggest that the upper mass cutoff may vary from 
galaxy to galaxy. 

Using different arguments based on various luminosity 
ratios, Deutsch & Willner 1986, Coziol & Demers 1995 and 
Coziol 1996 have found that B type stars may be predom- 
inant in the ionised regions of massive starburst galaxies. 
Taken at face value, this last phenomenon seems consistent 
with the ages of bursts deduced by Goldader et al. and the 
special form of the IMF found by these authors and Doyon, 
Puxley & Joseph. However, as pointed out by Coziol 1996, 
assuming starburst is a self-sustained phenomenon, this ob- 
servation could also be interpreted in terms of a sequence of 
short bursts spread over a few Gyr period. 

Evidence of extended or multiple bursts was recently 
found in many well known active star forming galaxies, like 
NGC 6764 (Schinnerer, Eckart & Boiler 2000), Arp 220 
(Anantharamaiah et al. 2000), NGC 1614 (Alonso-Herrero 
et al. 2000) and most recently M82 (de Grijs, O'Connell & 
Gallagher 2000). In another recent paper, Meurer (2000), 
using HST, distinguished in nearby starburst galaxies dif- 
ferent star clusters embedded in a diffuse glow of recently 
formed stars. Examining their colors, he concluded that 
these clusters are consistent with instantaneous bursts with 
very young ages (from up to 100 Myr), while the diffuse 
light seems to be produced by stars formed continuously 
over a period ten times longer than the crossing times. A re- 
cent study of clusters in the irregular starburst galaxy NGC 
1569 reveals similar clusters with ages between 3 Myrs and 
1 Gyr (Hunter et al. 2000). Old (3 Gyr) merger remnants 
were recently reported in the radio galaxy NGC 1316 (Goud- 
frooij et al. 2000). Multiple bursts and low upper mass cut- 
off for the IMF may also be necessary to explain the num- 
ber of different Wolf-Rayet sub-types in some WR galaxies 
(Schaerer et al. 2000). Even the small mass and less evolved 
HII galaxies seem to be formed of age-composite stellar sys- 
tems (Raimann et al. 2000). 

In the present article, we explore these issues further 
by studying the spectral characteristics in the near infrared 
(NIR) of a new sample of 48 starburst galaxies with differ- 
ent physical characteristics. As it is well known, the star- 
burst family is composed of a great variety of galaxy types, 
spanning several orders in size, mass, luminosity and gas 
metallicity (Salzer, MacAlpine, & Boroson 1989). The burst 
characteristics of these galaxies may not necessarily be the 
same. By selecting galaxies with different physical charac- 
teristics, we therefore expect to draw a more general picture 
of the starburst phenomenon. 

Our observations complete the study of luminous in- 
frared galaxies (LIRG) made by Goldader et al. (1995, 
1997a,b). By adding these two samples together, one can 
compare starburst galaxies which span a factor 100 in B 
and almost 1000 in FIR luminosity. 

All the galaxies in our sample have modest far infrared 
(FIR) luminosity (L IR < 10 10 L Q ). According to the IRAS 
luminosity function (Rieke & Lebofsky 1986; Vader & Simon 
1987; Saunders et al. 1990) such galaxies may form the bulk 
of the star forming galaxies in the nearby Universe. Our 
study, therefore, should reveal the "normal" behaviour in 
the NIR of starburst galaxies in the nearby Universe. 



The K-band window in the NIR is ideal for our study. 
The extinction by dust in this part of the spectrum is ~ 1/10 
that in the optical, which allows to compare star formation 
in galaxies with very different dust content. The K-band 
window also offers several diagnostic spectral lines for study- 
ing young stellar populations: the narrow emission line Br7 
provides some constraints on the relative number of O and 
B stars, while the CO band absorption, longward of 2.3 /im, 
can be used to estimate the ratio of Red Giants (RGs) over 
Red Super-Giants (RSGs) and better constrain the duration 
and age of the bursts (Doyon, Joseph & Wright 1994). All 
these parameters can now be easily deduce from available 
synthetic starburst models (Leitherer et al. 1999). 



2 PRESENTATION OF THE SAMPLE 

2.1 Selection criteria and biases 

The 48 galaxies in our sample were selected using optical 
criteria: they present a UV-excess in spectral lines or in the 
continuum, as established using an objective prism or the 
multiple filters technique. The galaxies were taken from var- 
ious sources: the Montreal Blue Galaxy (MBG) survey (Co- 
ziol et al. 1993; 1994), the sample of compact Kiso galaxies 
observed by Augarde et al. (1994), the University of Michi- 
gan survey (Salzer, MacAlpine, & Boroson 1989) and the 
spectrophotometric catalogue of HII galaxies (Terlevich et 
al. 1991). 

Discussion of the biases affecting these various surveys 
can be found in Coziol et al. (1997). In principle, by adding 
galaxies from objective prism and multiple filters surveys 
one can span the complete spectrum of starburst galaxies, 
including HII galaxies and Starburst Nucleus Galaxies (SB- 
NGs). In practice, however, observation of HII galaxies in 
the NIR is more difficult, being less luminous than the SB- 
NGs, the HII galaxies need too long exposure times. This 
explain the bias towards SBNGs in our sample. 

By choosing galaxies which were selected based on op- 
tical criteria, we also favour galaxies which have a moderate 
to low FIR luminosity (Coziol et al. 1997). As we already 
mentioned in our introduction, this bias is opportune, as it 
allows us to study the most common type of starburst galax- 
ies found in the nearby Universe (Rieke & Lebofsky 1986; 
Vader & Simon 1987; Saunders et al. 1990). 

2.2 Physical characteristics 

The physical characteristics of the galaxies in our sample are 
listed in Table 1. The redshifts, B magnitudes, morphologies 
and projected dimensions in the optical were all taken from 
NEDfj. The distances, in column 3, were estimated assuming 
Ho = 75 km s _1 Mpc" 1 , after correcting for the motion of 
the sun. These distances were used to determine the linear 
scale in kpc subtended in each galaxy by an angle of 1 arc- 
min on the sky (column 7). In many cases the morphology 
was not given. An examination of their images, as available 
in NED, allows us to determine that they look compact in 
appearance. These cases are identified in column 5. Only 

* The NASA/IPAC Extragalactic Database. 
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Table 1. Characteristics of the sample 
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one galaxy, MBG 23388-1514, was not morphologically clas- 
sified. Although we did not judged it compact, we could not 
establish its morphology. 

Various spectroscopic data were collected from the lit- 
erature to complement our analysis. Optical spectra of these 
galaxies were obtained using comparable apertures (2-3 
arc-second) and median-low spectral resolutions, allowing 
fair comparisons with our data. When the flux is avail- 
able, the observed Ha luminosity (column 8) was estimated. 



Metallicities (column 9) were determined using the ratio 
R3 = 1.35x([OlIl]A5007/H/3) and the empirical relation 
proposed by Vacca & Conti (1992). For the metallicities of 
I Zw 18 and II Zw 40, we adopted the values estimated by 
Izotov, Thuan & Lipovetsky (1998) and Masegosa, Moles 
& Campos-Aguilar (1994). The activity types (column 10) 
were established by comparing the ratios [Nil] A6584/Hq 
with [OIII]A5007/H/3 (Baldwin, Phillips & Terlevich 1981; 
Veilleux & Osterbrock 1987). The numbers in the last col- 
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Table 2. FIR characteristics for some galaxies in our sample 
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umn identify the references for the various spectroscopic 
data. 

We give in Table 2 some information in the FIR for 
31 galaxies (64%) in our sample. These data come from 
the IRAS Faint Source Catalogue, as found in NED. The 
FIR luminosity in column 2 was determined using the 
relation (Londsdale et al. 1985): log(Lm) = log(FiR) + 
21og[z(z + 1)] + 57.28, where z is the redshift, Fir = 
1.26 x 10- n (2.58/ 6 o + /too) erg cm" 2 s" 1 and / 60 and / 100 
are the fluxes in Jansky at 60 and 100 /im respectively. 
In column 3 and 4, we also give the FIR spectral indices 
(Sekiguchi 1987): a{\ x ,\-x) = log(fxJfx 2 )/ log(A 2 /Ai). To 
facilitate our analysis, only galaxies which have IRAS fluxes 
with high or intermediate qualities had their FIR character- 
istics estimated. The last column gives the classification of 
activity type as deduced from FIR criteria (see section 2.2). 
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2.3 Nature of the starburst galaxies in our sample 
and physical variety of starburst galaxy hosts 

We show in Figure la the spectroscopic diagnostic diagram 
used to classify 28 galaxies in our sample with spectroscopic 
data available in the optical. In this diagram the horizon- 
tal dash line distinguishes between HII galaxies and SBNGs 
(Coziol 1996). We count 19 SBNGs for 9 HII galaxies. This 
result confirms our observational bias towards SBNGs. The 
thin continuous curve in Figure la follows the locus traced 
by normal disc HII regions (Coziol et al. 1994). On this 
curve, the gas metallicity, given as 12 + log(0/H), increases 



Figure 1. a) Spectroscopic diagnostic diagram in the optical. 
The bold dash curve separates starbursts from AGNs. The hori- 
zontal dash line separates high from low excitation emission-line 
galaxies. The thin continuous curve is the HII regions sequence in 
metallicity (solar metallicity = 8.84 on this scale), b) IRAS diag- 
nostic diagram with two blackbodies model for dust. The numbers 
indicate the fractional contribution from hot dust. The horizon- 
tal dash line separates SBNGs from AGNs. The gray area is the 
locus occupied by HII galaxies, c) FIR vs B luminosities. The 4 
diagonals correspond to different ratios Lir/Lb 
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as the level of excitation (related to [OIIIJ/H/9) decreases. 
The galaxies in our sample cover the complete sequence in 
starburst's metallicities. Considering only the SBNGs, we 
found a mean metallicity 12 + log(0/H) = 8.7 ± 0.2, which 
is typical of this class of galaxies (Coziol et al. 1999). 

Half of the galaxies that could not be classified in the 
optical have FIR characteristics which can be used to deter- 
mine their starburst type: SBNGs or HII galaxies. In Fig- 
ure lb we show a diagnostic diagram based on IRAS spec- 
tral indices a(60,25) and a(100,60). Starburst galaxies are 
described in this diagram by a model with two dust com- 
ponents (Sekiguchi 1987): a cold component with a dust 
temperature of 27 K and a hot component with dust tem- 
peratures varying between 70 and 90 K. According to this 
model, starburst galaxies differ from normal spiral galaxies 
by their hot dust contribution, which is usually higher than 
60%. Although very few HII galaxies are detected in the 
FIR, they can usually be distinguished from SBNGs from 
their hotter dust temperatures. This is indicated by the gray 
area in Figure lb. According to this last criterion, most of 
the galaxies with FIR data in our sample are SBNGs. 



2.4 Probability to find an AGN in our sample 

Using FIR criteria we can also ascertain what is the prob- 
ability that one of the galaxies in our sample is dominated 
by an AGN. According to Coziol et al. (1998), 99% of IRAS 
starburst galaxies have a(60,25) < —1.9 (the horizontal dash 
line in Figure lb), while only 25% of Seyfert galaxies have 
such FIR colours. Almost all the galaxies in our sample have 
a(60,25) < —1.9, which excludes the possibility that one of 
these galaxies is dominated by an AGN. Note that the three 
galaxies exceeding this limit in our sample (Mrk 538, Mrk 
401 and MBG 02072-1022) have already been classified as 
SBNGs from their spectra. In fact, all the galaxies previously 
classified as starburst based on their optical spectra occupy 
this region of the diagram, which thus confirms their nature 
as starburst. 

Only six of the galaxies detected in the FIR diagram 
do not have a spectral classification in the optical. Based 
on the reported statistics above, and the total number of 
galaxies which was necessary to established it (see Coziol 
et al. 1998), we estimate to 7% the probability to find an 
AGN (which is probably not predominant) in this region of 
the diagram. This reduce to none the probability that one 
of these six galaxies also host an AGN. This results is not 
surprising considering the moderate FIR luminosity of these 
galaxies and the fact that the probability to find an AGN 
increases with the FIR luminosity (Veilleux, Kim & Sanders 
1999). 

The above conclusion hold also for the ten galaxies in 
our sample without spectral classification or FIR detection. 
Height of these galaxies were classified as starburst based on 
an objective prism spectrum, which rules out the presence 
of a broad line or luminous AGN. The last two are MBG 
galaxies which were classified as starburst based on narrow 
lines and small [Nil] /Ha ratios in one spectrum. The fact 
that these galaxies were not detected in IRAS suggests that 
the probability they host an AGN is practically null. 



2.5 Comparison with LIRGs and uLIRGs 

At this point of our analysis, it is convenient to compare 
our sample of galaxies with the 16 LIRGs previously classi- 
fied as starburst in the study of Goldader et al. (1997a). In 
Figure la, we see that the LIRGs have spectroscopical char- 
acteristics similar to SBNGs. However, they are generally 
more metal rich, having solar metallicity on average. Their 
hot dust contribution (Figure lb) is also higher than in SB- 
NGs, which is consistent with their higher FIR luminosity. 

In Figure lc, we compare the FIR and B luminosities 
of our galaxies with those of the LIRGs. For comparison, we 
also include some ultra luminous infrared galaxies (uLIRGs) 
from the literature. The ratios Lir/Lb can be used to distin- 
guish galaxies with unusually high star formation rates (see 
Coziol 1996): normal late-type spiral galaxies with constant 
star formation rates over the last 3 Gyrs have 1/3 <Lir/Lb 
< 3. One can see that the ratio Lir/Lb tends to increase 
with the B luminosity. The uLIRGs and LIRGs have both 
higher B luminosities and, therefore, higher Lir/Lb ratios 
as compared to SBNGs. Most of the galaxies in our sample 
have ratios Lir/Lb varying between 1 and 3. This compar- 
ison suggests that uLIRGs and LIRGs are generally more 
massive than SBNGs (having higher B luminosity) and may 
have more active burst regions or have younger bursts (as 
judged from their higher Lir/Lb ratios). 

One remarkable feature of Figure lc is how the starburst 
galaxies in our sample differ from well known uLIRGs, like 
Arp 220, NGC 6240, NGC 1614 or NGC 3690. The nature 
of uLIRGs is still in debate. It is not sure if part of the 
activity observed in these galaxies is not due to an hidden 
AGN. This is not the case of the galaxies in our sample. The 
position, in Figure lc, of the "prototype" starburst galaxy 
M82 seems, however, interesting. This galaxy looks more like 
an uLIRG than to one of our galaxies. This is consistent with 
the extreme properties observed in this object, which may 
also suggest some kind of AGN activity (Seaquist, Frayer 
& Frail 1997; Allen & Kronberg 1998; Kaaret et al. 2000). 
Two other archetype starburst galaxies in our sample are 
NGC 7714 (Mrk 538) and Mrk 710. These two galaxies are 
more representative of the starburst activity in our sample 
of galaxies. 

The above analysis emphasizes the diversity of physical 
characteristics of galaxies hosting a starburst. Our observa- 
tions should reveal if this variety of physical characteristics 
is reflected in their NIR spectral features. 



3 OBSERVATION 

The 2-2.5 /im spectra, studied in this paper, were obtained 
with the CGS4 spectrometer (Mountain et al. 1990) on 
UKIRT on the nights of 1992 November 5 and of 1993 
November 24, 25, and 26. A 75 1/mm grating was used in 
first order, providing a resolving power (X/5X) of ~300 in 
the K-ba,nd and a spectral coverage encompassing most of 
the if-band. The plate scale was 3 arc-second per pixel and 
the slit width was 1 pixel. This aperture was sufficient to 
cover the central and most luminous part (usually the nu- 
cleus) of the galaxies. The slit length was 30 pixels. For all 
galaxies, the slit was positioned in the north south direction 
on the brightest optical peak as seen from the offset guider 
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Table 3. Fluxes and equivalent widths for Br7 and H2 lines 



Aperture 1 Aperture 2 



Name 


z 


f(Br 7 ) 


EW(Br 7 ) 


f(H 2 ) 


EW(H 2 ) 


f(Br 7 ) 


EW(Br 7 ) 


f(H 2 ) 


EW(H 2 ) 






xl0~ 15 




xlO" 15 




xlO -15 




xlO -15 








(erg s _1 ) 


(A) 


(erg s" 1 ) 


(A) 


(erg s" 1 ) 


(A) 


(erg s" 1 ) 


(A) 


MBG00027-1645 


0.0269 


2.1 ± 0.3 


4.1 ± 0.7 


1.7 ± 0.3 


3.2 ± 0.6 


0.7 ± 0.2 


2.6 ± 0.8 


< 0.8 


< 2.9 


MBG00439-1342 


0.0019 


< 2.7 


< 3.7 


< 1.6 


< 2.1 


< 0.9 


< 2.7 


< 0.7 


< 2.1 


MBG00461-1259 


0.0232 


1.8 ± 0.4 


11.0 ± 2.7 


< 0.8 


< 4.5 


0.7 ± 0.2 


8.3 ± 1.7 


< 0.3 


< 3.2 


UM306 


0.0175 


< 0.2 


< 10.9 


< 0.2 


< 7.2 


< 0.1 


< 11.2 


< 0.1 


< 11.6 


Mrkl002 


0.0121 


8.8 ± 0.6 


8.4 ± 0.6 


2.7 ± 0.8 


2.5 ± 0.7 


5.5 ± 0.4 


8.6 ± 0.6 


< 1.3 


< 1.9 


UM372 


0.0395 


< 0.9 


< 11.7 


< 0.7 


< 8.4 


< 0.4 


< 5.4 


< 0.5 


< 6.7 


Mrk363 


0.0087 


< 1.0 


< 3.0 


< 1.6 


< 4.4 


< 0.4 


< 2.7 


< 0.6 


< 3.8 


MBG02072-1022 


0.0138 


37.4 ± 1.9 


10.5 ± 0.5 


< 11.5 


< 3.1 


25.2 ± 0.8 


13.9 ± 0.4 


< 5.4 


< 2.8 


MBG02141-1134 


0.0149 


4.2 ± 0.5 


7.7 ± 0.9 


2.0 ± 0.4 


3.5 ± 0.6 


1.7 ± 0.4 


7.6 ± 1.7 


0.8 ± 0.2 


3.5 ± 0.8 


Mrkl055 


0.0388 


< 1.1 


< 4.2 


< 1.1 


< 3.7 


< 0.8 


< 4.1 


< 0.6 


< 2.8 


MRK602 


0.0018 


< 6.1 


< 7.3 


< 3.0 


< 3.5 


< 5.5 


< 9.1 


< 2.2 


< 3.5 


KUG0305-009 


0.0214 


1.1 ± 0.3 


4.2 ± 1.3 


< 0.8 


< 2.8 


0.6 ± 0.2 


3.9 ± 1.2 


< 0.6 


< 3.8 


Mrk603 


0.0095 


23.5 ± 0.8 


22.9 ± 0.8 


4.1 ± 0.8 


3.7 ± 0.8 


16.9 ± 0.4 


23.4 ± 0.6 


1.4 ± 0.3 


1.8 ± 0.4 


KUV03073-0035 


0.0236 


1.0 ± 0.3 


5.0 ± 1.4 


< 1.0 


< 4.5 


0.6 ± 0.1 


5.6 ± 0.9 


< 0.6 


< 5.1 


MBG03084-1059 


0.0180 


< 1.2 


< 4.3 


< 0.9 


< 3.3 


0.8 ± 0.2 


5.4 ± 1.1 


< 0.5 


< 3.3 


MBG03183-1853 


0.0132 


< 0.7 


< 3.0 


< 1.1 


< 4.4 


< 0.4 


< 3.5 


< 0.4 


< 3.8 


MBG03317-2027 


0.0047 


< 1.6 


< 2.1 


< 2.8 


< 3.5 


< 1.4 


< 2.8 


< 1.3 


< 2.4 


KUG0338+032 


0.0228 


< 0.7 


< 8.3 


< 0.8 


< 8.2 


0.4 ± 0.1 


5.8 ± 1.8 


<0.3 


< 4.4 


MBG03468-2217 


0.0157 


8.6 ± 0.7 


8.4 ± 0.7 


3.3 ± 0.5 


3.1 ± 0.4 


6.3 ± 0.5 


8.2 ± 0.6 


1.6 ± 0.4 


2.0 ± 0.5 


MBG03523-2034 


0.0070 


1.3 ± 0.4 


5.4 ± 1.7 


1.1 ± 0.3 


4.7 ± 1.1 


< 0.5 


< 3.7 


< 0.3 


< 2.0 


IRAS04493-0553 


0.0095 


1.6 ± 0.3 


7.4 ± 1.4 


0.9 ± 0.3 


3.9 ± 1.2 


0.5 ± 0.2 


6.9 ± 2.3 


< 0.4 


< 5.5 


Mrkl089 


0.0148 


10.7 ± 0.3 


84.9 ± 3.9 


1.8 ± 0.2 


13.5 ± 1.9 


3.7 ± 0.2 


79.0 ± 4.7 


0.9 ± 0.2 


18.8 ± 4.5 


Mrkll94 


0.0167 


5.6 ± 0.5 


4.1 ± 0.4 


3.7 ± 0.6 


2.6 ± 0.4 


1.1 ± 0.3 


1.8 ± 0.5 


1.5 ± 0.3 


2.3 ± 0.4 


IIZW40 


0.0034 


45.1 ± 0.4 


264.4 ± 9.3 


2.3 ± 0.5 


13.0 ± 2.7 


34.3 ± 0.4 


322.0 ± 15.0 


1.2 ± 0.2 


11.6 ± 1.4 


KUG0720+335 


0.0134 


16.9 ± 0.6 


18.4 ± 0.7 


3.5 ± 0.7 


3.6 ± 0.8 


13.2 ± 0.4 


18.9 ± 0.6 


2.1 ± 0.4 


2.8 ± 0.5 


Mrk384 


0.0155 


3.0 ± 1.0 


3.7 ± 1.2 


< 3.1 


< 3.7 


2.0 ± 0.4 


5.0 ± 0.9 


< 1.5 


< 3.5 


KUG0815+249 


0.0075 


1.5 ± 0.4 


7.7 ± 1.9 


< 1.0 


< 4.8 


0.6 ± 0.1 


5.8 ± 1.1 


< 0.7 


< 6.2 


KUG0815+246 


0.0091 


< 0.9 


< 7.4 


< 1.4 


< 10.4 


< 0.6 


< 8.5 


< 0.7 


< 8.7 


KUG0821+229 


0.0253 


0.5 ± 0.2 


4.2 ± 1.4 


< 0.9 


< 7.1 


0.5 ± 0.2 


8.1 ± 2.5 


< 0.8 


< 11.1 


KUG0825+252 


0.0071 


< 0.4 


< 3.9 


< 0.7 


< 6.5 


< 0.3 


< 5.5 


< 0.3 


< 4.7 


Mrk90 


0.0112 


< 2.0 


< 4.8 


< 1.0 


< 2.4 


< 0.8 


< 4.2 


< 1.0 


< 4.7 


Mrkl02 


0.0105 


< 2.0 


< 3.5 


< 2.5 


< 4.1 


1.2 ± 0.4 


4.4 ± 1.5 


< 1.8 


< 6.2 


Mrk401 


0.0061 


6.8 ± 0.7 


5.2 ± 0.5 


< 2.4 


< 1.7 


6.4 ± 0.5 


7.3 ± 0.6 


< 1.5 


< 1.6 


IZW18 


0.0020 


< 2.0 


< 180.4 


< 1.3 


< 85.3 


< 0.8 


< 76.5 


< 0.7 


< 70.4 


Mrk404 


0.0038 


11.5 ± 1.3 


5.3 ± 0.6 


< 4.4 


< 1.9 


8.9 ± 0.7 


6.1 ± 0.5 


3.0 ± 1.0 


2.0 ± 0.6 


Mrk710 


0.0055 


18.3 ± 0.6 


36.3 ± 1.2 


3.1 ± 0.6 


5.7 ± 1.1 


15.1 ± 0.4 


57.5 ± 2.2 


2.3 ± 0.3 


8.3 ± 1.2 


Mrk33 


0.0040 


29.8 ± 0.9 


34.0 ± 1.1 


3.2 ± 0.8 


3.5 ± 0.9 


18.6 ± 0.5 


33.8 ± 1.0 


2.0 ± 0.4 


3.5 ± 0.7 


MBG21513-1623 


0.0381 


0.5 ± 0.1 


6.5 ± 1.5 


< 0.6 


< 7.0 


< 0.4 


< 11.3 


< 0.3 


< 7.1 


Mrk307 


0.0193 


< 2.1 


< 8.9 


< 1.4 


< 5.4 


1.2 ± 0.3 


10.4 ± 2.3 


< 0.8 


< 6.3 


MBG22537-1650 


0.0138 


1.0 ± 0.3 


3.5 ± 1.0 


< 1.2 


< 4.0 


0.6 ± 0.2 


4.0 ± 1.0 


< 0.6 


< 3.5 


KUG2254+124 


0.0256 


0.8 ± 0.2 


10.6 ± 2.7 


< 0.8 


< 9.9 


0.4 ± 0.1 


8.8 ± 2.8 


< 0.5 


< 10.8 


KUG2300+163 


0.0082 


1.3 ± 0.2 


21.4 ± 2.5 


< 0.5 


< 7.0 


0.6 ± 0.1 


22.0 ± 2.8 


< 0.2 


< 5.9 


Mrk326 


0.0140 


5.2 ± 0.6 


5.3 ± 0.6 


< 2.4 


< 2.3 


6.1 ± 0.4 


8.4 ± 0.6 


2.4 ± 0.6 


3.1 ± 0.8 


MBG23318-1156 


0.0206 


0.9 ± 0.2 


12.9 ± 3.0 


< 1.0 


< 12.4 


0.5 ± 0.2 


15.4 ± 4.9 


< 0.4 


< 10.1 


Mrk538 


0.0110 


14.6 ± 0.9 


17.1 ± 1.1 


6.7 ± 1.0 


7.5 ± 1.1 


10.5 ± 0.4 


18.9 ± 0.7 


3.2 ± 0.3 


5.4 ± 0.6 


MBG23372-1205 


0.0239 


0.8 ± 0.2 


9.0 ± 2.7 


< 0.7 


< 7.6 


< 0.9 


< 21.3 


< 0.8 


< 16.8 


MBG23388-1514 


0.0303 


< 1.6 


< 23.3 


< 1.6 


< 23.0 


< 0.5 


< 13.7 


< 0.6 


< 14.8 


Mrk332 


0.0110 


2.5 ± 0.8 


2.6 ± 0.8 


2.6 ± 0.7 


2.6 ± 0.7 


2.2 ± 0.6 


3.7 ± 0.9 


1.9 ± 0.5 


2.9 ± 0.8 



camera. Integration times were generally 10-20 minutes and 
sky frames were obtained by nodding 10 pixels along the slit, 
hence keeping the object always on the slit. Several spectro- 
scopic and photometric standards were also observed during 
each night. The spectrophotometric accuracy of our spectra 
is estimated to be ~20%. The spectra were reduced under 
the IRAF environment using IRAF scripts written by one of 
us (RD). Details on the procedure followed for the reduction 
can be found in Doyon et al. (1995). 



In Table 3, we present the fluxes and equivalent widths 
of the Br 7 and H 2 emission lines as measured in our spectra. 
These quantities were estimated after rectifying the contin- 
uum by fitting a power law (F\ oc A* 3 ) from the featureless 
sections of the spectrum. The choice of a power-law is physi- 
cally justified by the fact that late-stars show approximately 
Rayleigh-Jeans spectrum in the NIR (Fa oc A~ 4 , between 
2 and 2.9/wn; see appendix A of Doyon, Joseph & Wright 
1994). A x 2 minimisation routine was then used to fit a 
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Table 4. Continuum slope, CO index and K magnitudes 
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Gaussian profile to the Br7 emission line. In Table 3, val- 
ues without an error estimate corresponds to 3 sigma upper 
limits. 

We give in Table 4 the exponent (/3) of the power law fit 
for the continuum, the CO spectroscopic indices (CO spe c), 
as defined in Doyon, Joseph & Wright (1994), and the spec- 
troscopic K magnitudes. The CO spec quoted include the un- 
certainty in the continuum level. In general, the error on 
CO spec is dominated by this uncertainty, which is less than 
0.02 mag. The second column of Table 4 gives the redshifts 
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as deduced from the features observed. These redshifts are in 
relatively good agreement with the one found in NED. The 
last column gives the signal to noise ratios in our spectra. 

For each galaxy in the two tables we give two sets of 
data. The first set corresponds to values as measured using 
an aperture covering 3x9 arc-second (Aperture 1). The sec- 
ond set (the last four columns of each tables) corresponds to 
values as measured using an aperture covering only 3x3 arc- 
second (Aperture 2). In our following analysis, the smaller 
aperture will be associated to the nucleus of the galaxies. 
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2.1 2.2 2.3 2.4 2.5 2.0 2.1 2.2 2.3 2.4 




2.1 2.2 2.3 2.4 2.5 2.0 2.1 2.2 2.3 2.4 




Rest wavelength (|im) Rest wavelength (|im) 

Figure 2. Some examples of spectra with strong Bracket 7 emission— line. The spectra were shifted to their rest frame using velocities 
in Table 1. Conspicuous emission lines and absorption features are identified. For each galaxy we show the power law fitted on the 
continuum and the window used to measure CO S pec- We also show a straight line fit. A straight line do not give a good fit for the K 
continuum of late-stars and galaxies. However, it would have no significant effects on the measurements of the CO spe c (variation < 1%) 
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Rest wavelength (|im) Rest wavelength (|im) 

Figure 3. Examples of spectra with relatively weak Bracket 7 emission-line and strong CO lines. The spectra are as described in Figure 2 



© 2000 RAS, MNRAS 000, 



10 R. Coziol, R. Doyon & S. Demers 



A few representative spectra are shown in Figure 2 and 
Figure 3. We regroup the spectra in two categories. In the 
first category, Figure 2, the galaxies possess strong Bracket 7 
(2.166 (im) emission. These galaxies also show relatively 
strong Hel (2.059 /Um) and H 2 1-0 S(l) (2.122 /im) emis- 
sion lines. They show, on the other hand, different strengths 
of the CO absorption band longward of 2.293 fim. 

The most extreme case in our sample is the HII galaxy 
II Zw 40. This galaxy shows very strong emission lines and 
extremely weak CO band. This galaxy also happen to be one 
of the less chemically evolved galaxies in our sample. Two 
SBNGs with intense emission lines are Mrk 710 and Mrk 
1089. Most galaxies in our sample resemble Kug 0720+335 
and Mrk 538 (NGC 7714), which have medium intensity 
emission lines and relatively strong CO band. 

In the second category, Figure 3, we find galaxies which 
have weaker emission lines than Mrk 538 and somewhat 
stronger CO band. In some galaxies, like Mrk 332, the NIR 
emission lines are barely visible, and only upper limits can 
be measured. 

The NIR spectra of our sample of galaxies show sig- 
nificant variations. This result indicates that the diversity 
of physical characteristics of starburst galaxy hosts, as ob- 
served in the optical and far infrared, implies some variety 
of the NIR spectral features. Synthetic starburst models will 
now be used in order to understand what these variations 
mean in terms of the burst characteristics. 



4 NIR ANALYSIS 

4.1 Synthetic starburst models: duration, age of 
the burst and IMF 

Using synthetic starburst models, it is possible to interpret 
NIR features like EW(Br7) and CO spcc in terms of the dura- 
tion of the burst, its age and the form of the IMF (Puxley, 
Doyon & Ward 1997; Leitherer & Heckman 1995; Mayya 
1997). The principal steps in the synthetic starburst model 
are the following (Leitherer & Heckman 1995): A) stars are 
formed at specified rates and distributed along the HR dia- 
gram; B) evolutionary models describe the time dependence 
of the individual physical parameters of the stars (mass, lu- 
minosity, etc.); C) the stellar number densities in the HR 
diagrams are determined and any desired synthetic quan- 
tities assigned to each stars. By summing over the entire 
stellar population the model yields the integrated properties 
of starburst galaxy. 

In these models, star formation is represented by an ex- 
ponential law: £(t) = £(i = 0) exp — t/r, where t is the age 
of the burst and r its duration. Two limiting cases are of 
special interest: an instantaneous burst and a constant star 
formation rate. An instantaneous burst corresponds to the 
case where the duration of star formation is short as com- 
pared to the age of the galaxy (r << t). All the stars are 
thus formed at the same epoch. The other limiting case is a 
constant star formation rate (CSFR). This mode of star for- 
mation is characteristic today of disc in giant spiral galaxies 
(Kennicutt 1983). 

It is important to note that theoretically a CSFR can 
always be approximated as a series of instantaneous bursts. 
It means that a system having experienced a sequence of 



bursts over a time period comparable to r = t would be 
observationally indistinguishable from a galaxy with a CSFR 
(Leitherer et. all 1999). 

Another important parameter in synthetic starburst 
model is the IMF. This function is usually expressed as a 
power law: 4>{M) — Cm~ a , where the constant C is de- 
termined by the total mass converted into stars. To get the 
whole spectrum of star masses formed, the IMF is integrated 
between the upper and lower mass cutoff, M up and M; ou) . In 
starburst models the exponent of the IMF and the two mass 
cutoff are not very well constrained and it is therefore im- 
portant to try different possible values. 

We choose to compare our data with the models of Lei- 
therer et al. (1999). These authors have made their grids of 
results directly accessible on the WEbQ We consider only 
two limiting cases: the instantaneous burst and the CSFR (1 
Mq per yr). For each of these cases, three different scenarios 
are tested: 1) a Salpeter IMF, a = 2.5, with M up = 100 M ; 
2) a Salpeter's IMF with M up = 30 M ; 3) an IMF with a 
steeper slope, a = 3.3 and M„ p = 100 Mq. The models 
cover ages from 10 6 to 10 9 yr. 

Comparisons of our data for EW(Br7) and CO spec with 
the instantaneous burst models are presented in Figure 4. 
From this figure, we conclude that none of the instanta- 
neous burst models can fit the ensemble of our data. For 
solar metallicity (Z = 0.02) the models generally underes- 
timate the value of EW(Br7). Some galaxies can be fitted 
by instantaneous bursts, but the predicted metallicity is 2 
times solar, which is much higher than the gas metallicity we 
have estimated in these galaxies. Our results confirm those of 
Goldader et al. (1997b), extending their conclusion to star- 
burst galaxies with different physical properties and FIR 
luminosity. It implies that the instantaneous burst scenario 
is not a good representation of the starburst phenomenon. 

Our data are compared with the CSFR models in Fig- 
ure 5. One can see in Figure 5a that the Salpeter's IMF 
model with M up = 100 Mq cannot fit the ensemble of our 
data. But a relatively good fit is obtained in Figure 5b, when 
the upper mass cutoff is equal to M up = 30 Mq . We see no 
improvement in Figure 5c in adopting an IMF slope a — 3.3. 

Once again our results confirm those of Goldader et al. 
(1997b), extending their conclusion to all the starburst types 
known. We conclude that the best parameters which de- 
scribe starburst galaxies in the nearby Universe are a CSFR 
and a Salpeter's IMF with an upper mass cutoff M up = 30 

Mq. 

According to this solution, the bursts have ages be- 
tween 10 Myr and 1 Gyr, the majority being older than 
100 Myr. Comparing with LIRGs, we found, on average, 
younger bursts for these galaxies, with ages lower than 100 
Myr. This is consistent with the difference in luminosity ra- 
tios noted in Section 2.5. The model reproduce correctly also 
the fact that LIRG are more metal rich than the SBNGs. 
The consistency with the metallicities deduced from their 
optical spectra (Section 2.2) is remarkably good. 

According to the model, HII galaxies must have younger 
bursts than the SBNGs. The metallicities predicted by the 
models are also in good agreement with those measured. 
Young bursts and low metallicities are also consistent with 

t http:/ /www. stsci.edu/science/starburst99 
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Figure 4. Comparison of our data with instantaneous burst mod- 
els, a) IMF with slope a = 2.5 and M up = 100 M ; b) a = 2.5 
and M up = 30 M ; c) a = 3.3 and M up = 100 M . The models 
cover ages from 10 6 to 10 9 yr. Open circle = SBNG, filled square 
= HII galaxy, star = LIRG 



the higher dust temperatures suggested by the FIR model 
in Section 2.3. 

In our sample, the galaxy which have the youngest burst 
is II Zw 40. The position of this galaxy in Figure 5a suggests 
that it may be possible to fit an IMF with upper mass cutoff 
M up = 100 M (although the predicted metallicity would 
then be wrong). The same conclusion may also apply to 
Mrk 1089. These observations may suggest varying IMF in 
some galaxies (Doyon, Puxley & Joseph 1992). 

From the above discussion, we conclude that the solu- 
tion of the synthetic starburst model is consistent with ob- 
servations in the optical and FIR. In particular, it satisfacto- 
rily explains the differences between the different starburst 
types: HII galaxies have younger bursts and lower metal- 
licities than SBNGs, while LIRGs have younger bursts but 
higher metallicities. These differences may be related to dif- 
ferent histories of star formation or different stages in burst 
evolution. 




0.4 0.8 1.2 1.6 2.0 2.4 2.8 
Log EW(Bry) (A) 



Figure 5. Comparison of our data with constant star formation 
models (1 M Q per yr). a) IMF with slope a = 2.5 and M up = 100 
M Q ; b) a = 2.5 and M up = 30 M ; c) a = 3.3 and M up = 100 
M . The signification of the symbols are the same as in Figure 4 

Taken at face value, the above result confirms the stud- 
ies of Doyon, Puxley & Joseph (1992), Goldader et al. 
(1997b) and Coziol (1996). The generality of our observation 
further suggests that this must be a characteristic of star- 
burst galaxies in the nearby Universe. However, the above 
solution may not represent strong constraints on the dura- 
tion of the burst or the IMF. This is because the synthetic 
starburst models we used cannot distinguished between real 
constant star formation and a series of instantaneous bursts 
(Leitherer et. all 1999). In order to verify this alternative, 
we will now reexamine our spectra and search for traces of 
previous bursts. 

4.2 Spatial variation of NIR features: evidences 
for sequences of bursts 

In order to reproduce CSFR using a sequence of bursts, one 
must assume that this sequence spreads over a relatively long 
period of time. The solution of synthetic starburst model 
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Figure 6. Comparison of the ratios of Br7 fluxes (a) and of of 
the equivalent widths (b), as measured in and out of the nu- 
cleus and as function of the redshift. In this analysis we associate 
the smaller aperture (3x3 arc-second) to the nucleus. No rela- 
tion with redshift is observed, which indicates that the increase 
of scale with distance, implied by our two fixed apertures, does 
not introduce any systematic variations. The galaxy showing the 
higher variation is Mrk 1194. No particular reasons explain this 
variation 



suggests a period extending over 10 8 up to 1 Gyr. If this 
sequence of burst is related to some self-regulated mecha- 
nism or to propagation of star formation, we may also as- 
sume that the different bursts do not happen in the same 
regions of space. This implies that the bursts regions should 
be heterogeneous in space and in time. With sufficient spa- 
tial resolution, we should then observe a variation in burst 
characteristics as we survey different regions. As a first test, 
therefore, we can check if a variation of aperture in our spec- 
tra induces some variation of the NIR spectral features. 

Before doing the aperture test, it is important to ver- 
ify that the obvious increase of scale with distance, implied 
by our two fixed apertures, do not introduce any system- 
atic variations, which could be misinterpreted as variation of 
burst characteristics. Associating the smaller aperture (3x3 
arc-second) to the nucleus, we verify in Figure 6 that the 
ratios of the Br7 fluxes and equivalent widths, as measured 
in and out of the nucleus, do not show any relation with the 
redshift. Similar negative results were also obtained for H2, 
the K magnitudes, the CO spec and for the slope of the con- 
tinuum {f3). These negative results suggest that what ever 
variations we may observed for these ratios, their origin must 
be intrinsic to each galaxy. 

Having verified that the two fixed apertures do not in- 
troduce spurious effects, we now proceed to examine how the 
different ratios varied in each galaxy. In Figure 7a, we first 
compare the ratios of the fluxes in Br7 with the ratios of 
the equivalent widths. In general, the ratios f(Br7) /f(Br7) jv 



are higher than one, which indicate that the bursts of star 
formation generally extend over a few kpc around the nu- 
cleus. The equivalent widths, on the other hand, are gen- 
erally wider in the nucleus, except in galaxies where the 
ratio f(Br7) /f(Br7) jv > 2, which correspond to extended 
bursts, where the equivalent width stays constant or in- 
creases out of the nucleus. According to synthetic starburst 
models, EW(Br7) is an indicator of the age of the burst. 
From the above observation we then deduce that in the 
galaxies where the star formation is more concentrated, the 
bursts are younger in the nucleus. The contrary seems to 
be true in galaxies where the bursts are more extended. No 
difference is observed between HII galaxies and SBNGs. 

The case for H2, Figure 7b, is slightly different. Very few 
galaxies in our sample show significant H2 emission. In these 
galaxies we find the same tendency than for Br7, which is the 
H2 emission generally spreads out of the nucleus. But the H2 
equivalent width tends to be wider out of the nucleus. In Fig- 
ure 7c one can see that in galaxies where the star formation 
is concentrated in the nucleus (f(Br7)/f(Br7)jv < 2), the H2 
emission is more intense than Br7 out of the nucleus. Only 
two galaxies, Mrk 1089 and 1194, show the reverse tendency. 
For the galaxies in our sample with significant amount of 
H2 emission, the H2 emission seems, therefore, to be always 
more intense and wide spread in regions where Br7 is less 
intense and less abundant. We verified also that 50% of the 
galaxies without H2 detection have detectable Br7 emission, 
which is consistent with the general tendency observed. The 
signification of this variation is not obvious, however, and its 
possible interpretation will therefore be kept for later during 
our analysis. 

For the K-band magnitude, we naturally expect it to 
increase as the number of stars included in the aperture in- 
creases. However, the mean variation observed in Figure 7d 
seem smaller than what we expected. While the gain in sur- 
face represented by the difference in aperture is a factor 3, we 
find a meager AK = —0.7, which corresponds to an increase 
in flux by only a factor 1.8. The K-band emission, there- 
fore, seems mostly concentrated in the nucleus. Indeed, one 
can see in Figure 7d that the difference in K-band magni- 
tude is smaller in galaxies where the bursts are concentrated 
(f(Br7)/f(Br7)jv < 2) than in galaxies where the bursts are 
extended. The K-band emission seems, therefore, to tightly 
follow the region of the bursts. This observation suggests 
that the stars emitting in the K band must be tightly re- 
lated to the bursts. 

No trend was observed for CO spe c or f3. These two pa- 
rameters show random variations with standard deviations 
±0.04 for the difference in CO sp ec and ±0.15 for the ratios 
of p. 

To illustrate what effects the variations of NIR features 
have on the burst characteristics of each galaxy, we show in 
Figure 8 the changes introduced by an increase in aperture 
in the burst model of Figure 5b. In this figure, an increase 
of the equivalent width transforms into a younger age while 
an increase of the CO spe c transforms into an increase of 
metallicity and a younger age. In Figure 8, one can see that 
some variations may imply significant differences in terms 
of the burst characteristics. No trend is observed, however, 
which suggests that the cause of these variations may depend 
on the particular history of star formation of each galaxy. 



© 2000 RAS, MNRAS 000, 



Near infrared spectroscopy of starburst galaxies 13 



DQ 
LJJ 

lr 

DQ 
LU 



2.0 



SBNGs 
HII galaxies 



o 



1.0 



o 



2.0 




LU 



LJJ 



1.0 



O 



o 



o' 

-0* 



o 



o 
o o 



1.0 3.0 5.0 

f(Bry)/f(Bry) N 



1.0 



f(H 2 )/f(H 2 ) 



5.0 



2/N 




-0.2 
-0.4 



— -0.6 
* -0.8 



-1.0 



-1.2 



- 




o 



o 



1.0 3.0 

f(Bry)/f(Bry) N 



5.0 



0.0 2.0 4.0 

f(Bry)/f(Bry) N 



Figure 7. Variation of the NIR features as measured in two different apertures. We compare the values measured in a 3x3 arc-second 
aperture (nucleus) with those in a larger 3x9 arc— second aperture. The vertical dashed line in a, b and c is the limit adopted for spatial 
resolution 
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Figure 8. Variation of the burst characteristics with a variation 
in aperture. The parameters of the model are the same as in 
Figure 5b 
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Higher spatial resolution will be necessary to confirm this 
result. 

From the above analysis we conclude that the variation 
of the NIR spectral features with aperture are consistent 
with a variation of burst characteristics in space. This re- 
sult supports our hypothesis that the bursts regions are not 
homogeneous in space. However, the absence of systematic 
trends and the lack of spatial resolution impede us to con- 
clude on the real cause of this variation and we must there- 
fore search for other evidence in favour of the multiple burst 
hypothesis. 



4.3 Luminosity— luminosity relations 

Our previous analysis suggests that the stars emitting in 
the K-band are tightly related to the burst population. We 
should then expect to find a good correlation between the K- 
band absolute magnitude and parameters link to the burst. 
This is important to verify, because the stars producing the 
K-band luminosity are probably the same stars which pro- 
duce the CO band. A correlation of K-band luminosity with 
parameters directly linked to the burst would thus favour 
massive and young stars, namely RSGs, over older stars like 
RGs. 

We see in Figure 9a that L(Br7) increases with the K- 
band absolute magnitude. A linear regression performed on 
our data (continuous curve in Figure 9a) gives a weak but 
statistically significant coefficient of correlation R= 67%. 
(with 33 points, the probability to obtain higher correlation 



Figure 9. Relation between Br7 and absolute K magnitude. 
The continuous curves are linear regressions fitted using only our 
galaxies. The dash curves are linear regression fitted adding the 
LIRGs to our galaxies. The coefficients of correlation are given in 
the text 



by chance is P 33 ,o.67 < 0.05%; Pugh & Winslow 1966). The 
correlation improves to 82% when we include the LIRGs 
(dash curve in Figure 9a). The properties of our galaxies are 
in good continuity with those of the LIRGs, which form the 
upper limit of the distribution. 

The correlation of L(Br7) with the K-band absolute 
magnitude confirms that the stars producing this luminos- 
ity are related to the burst population. The relative weak- 
ness of this relation, however, may suggest that they are not 
correlated to the younger stars. But this result is consistent 
with what we expect if the stars emitting in the K-band are 
in majority RSGs: the number of RSGs increases with the 
intensity of the burst, but these stars also appear in great 
number only after the peak intensity of the bursts. 

In Figure 9b, we find a reverse tendency for EW(Br7), 
this parameter decreasing as increases. These two pa- 
rameters are not strongly correlated (R= 57%), although the 
correlation is still statistically significant (P33,o.57 < 0.05%). 
The obvious interpretation for this reverse tendency is that 
galaxies which are richer in burst regions naturally also con- 
tain the highest number of evolved stars. 

What is interesting in this last figure is the high dis- 
persion of the data. This high dispersion suggests that the 
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Figure 10. Relation between H2 and K-band absolute magni- 
tude. The signification of the symbols and curves are as in Fig- 
ure 9 



evolution of the bursts depends on the particular star for- 
mation history of each galaxy. This is consistent with our 
analysis of the variation with aperture and suggests that 
the bursts may also be heterogeneous in time (the disper- 
sion in EW(Br7) is a dispersion in time). The behaviour 
of the LIRGs is particularly remarkable, suggesting younger 
bursts than the observed tendency would have predicted. 
This could be one indication of repetitive bursts. 

Note that we do not reject, a priori, the alternative 
which is that the stars emitting in the K-band are old RGs. 
In this case, the correlation of Mk with L(Br7) would sug- 
gest some relation between mass and the intensity of the 
bursts, while the correlation between Mk and EW(Br7) 
would mean that the bursts evolve more rapidly in galaxies 
where they are more intense. However, it would be difficult in 
this case to understand why the correlation are not stronger. 
But such interpretations, in fact, are much more constrain- 
ing than what our data can revealed (we do not have infor- 
mation on the masses or the intensity of the bursts, which, 
anyway, should be normalized by the surface or volume). To 
our knowledge, such relations were never reported before, 
and our data are obviously insufficient to defend them. Our 
favoured interpretation, on the other hand, do not imply any 
new mechanisms, relying only on standard star evolutionary 
scenarios. Another argument in favour of our interpretation 
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Figure 11. Relations between Br7 and H2 luminosities. The sig- 
nification of the symbols are as in Figure 9. The diagonals corre- 
spond to ratios 1.0, 0.5 and 0.1, consistent with values for shocks 
or fluorescence. The short continuous curve is a linear regression 
on our galaxies. The dash curve is a linear regression including 
LIRGs 



may be the relatively strong correlation we find between 
L(H"2) and Mk, as we now show. 

A very strong correlation, R= 92%, is found between 
L(H2) and Mk in Figure 10a. The lower number of points 
in this fit do not influence the results: with 16 points, 
the probability to obtain a higher correlation by chance is 
Pi6,o.92 < 0.05%. This good correlation do not depends on 
the inclusion of the LIRGs either, the coefficient increasing 
only to R= 94% when we include them. For the equivalent 
widths, Figure 10b, we also find a better correlation than for 
Br7, R= 62%, Pi6,o.62 < 1.4%. This behaviour suggests that 
the H2 emission is better correlated to a more evolved phase 
of the bursts, being correlated to SRGs, than to a younger 
phase. In general, therefore, the behavior of H2 emission 
seems consistent with the hypothesis of heterogeneous burst 
in space and time, showing variation with aperture (section 
4.2) and the coexistence of evolved components with younger 
ones. It would help, however, to determine what is the pro- 
cess producing this line, since the two most probable mech- 
anisms, supernovae shocks remnants and UV-fluorescence, 
imply much different time scales. 

In principle, it would be possible to determine what is 
the mechanism producing the H2 emission by comparing the 
two luminosities L(Br7) and L(H2). In Figure 11 we find that 
L(H2) is relatively well correlated to L(Br7) (R= 74% with 
Pi6,o.62 < 0.4%). But the ratios L(Br7)/L(H2) stay between 
0.1 and 1, which is consistent with supernovae shocks rem- 
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4.4 Spectral types of ionising stars 
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Figure 12. a) Relation between Ho and Br7 luminosities. The 
different diagonals correspond to theoretical ratios between the 
luminosities, assuming various magnitudes of extinction as in- 
dicated by the numbers; b) Relation between the FIR and Br7 
luminosities. The diagonals are theoretical ratios expected for dif- 
ferent spectra types. The luminosities of Br7 were corrected for 
extinction as determined in a 



nants (Moorwood & Oliva 1990) or UV-fluorescence (Pux- 
ley, Hawarden & Mountain 1990). Therefore, although our 
observations suggest that the H2 emission phenomenon is 
related to a more evolved phase of the bursts, we are un- 
able to say if this phase happen a few million years (the 
time scale of supernovae) or a few 10 7 or 10 8 yrs (the time 
scale of B stars, responsible for UV-fluorescence) after the 
maximum of the bursts. 



In Coziol & Demers 1995 and in Coziol (1996), one evidence 
in favor of a sequence of bursts was the predominance of B 
type stars in starburst nucleus galaxies. We can repeat the 
analysis performed in Coziol (1996) to verify this result in 
the NIR. The advantage of working in the NIR is that the 
extinction in this part of the spectrum is less severe, and the 
fraction of obscured stars should therefore be less important 
(Calzetti et al. 1995), decreasing the uncertainties related to 
this method. 

The level of extinction in our galaxies can be estimated 
(Figure 12a) by comparing L(Br7) with the luminosity in 
Ha (L(Ha)). For the theoretical relation between L(Br7) 
and L(H«) we use the one adopted by Leitherer & Heck- 
man (1995). The theoretical ratio is traced in Figure 12a 
as a continuous curve. Lower ratios then must correspond 
to various magnitudes in extinction in the optical. In our 
galaxies, the extinction varied between and 2 magnitudes 
, which is in good agreement with values estimated from the 
Balmer decrement method (Contini, Considere & Davoust 
1998). This should not be surprising considering the low 
dust extinction in these galaxies. The LIRGs seem to suffer 
slightly higher extinction, the values varying between 3 and 
4 magnitudes. 

Once corrected for extinction, the flux ratios Br7/LiR 
are proportional to the ionising flux produced by stars with 
various spectral types (Devereux & Young 1990). We show in 
Figure 12b the locus expected for ionising star clusters with 
different spectral types. This figure confirms the dominance 
of B type stars as found in Coziol (1996). 



5 DISCUSSION 

The question that we need to answer is: can we accept the 
result of the synthetic starburst model? In order to answer 
this question we need first to better understand this result. 
The main reason why the instantaneous burst solution is re- 
jected by our observation is because the predictions of the 
models are in contradiction with the values of the CO spec- 
tral indices we measured. In an instantaneous burst, the 
Br7 equivalent width stays high for the first 10 7 years. Dur- 
ing this time, CO spec is negligible and begins to increase 
only around 10 7 year, when the most massive stars evolve 
as RSGs. It is because we observe both intense Br7 emis- 
sion and relatively strong CO lines that we need continuous 
star formation. The answer seems, therefore, to rely on the 
validity of the CO spectral index. 

To obtain a good fit with any of the instantaneous burst 
models, the CO spectral indices would need to decrease by 
more than 50%. Obviously, this cannot be due to the method 
we used to measure these lines. The strength of the CO 
bands in our spectra is undeniable. A relatively strong value 
of CO S pec is consistent with a predominance of RSGs over 
RGs (Doyon, Joseph & Wright 1994). This is consistent with 
the correlation found between Mk and L(Br/y), which sug- 
gests that the stars producing this luminosity, and probably 
also the CO band, are related to the burst population. More- 
over, this relation is relatively weak because RSGs appear 
in great number only slightly after the culmination of the 
burst. The relatively strong CO bands in our spectra seems 
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therefore to be real, and a constant star formation rate is 
consequently necessary to explain our data. 

But what about the low upper mass cutoff? If we look 
at Figure 5, we see that the reason why we need such feature 
is because the Br7 equivalent widths seem too weak as com- 
pared to the model with upper mass cutoff M up = 100 M Q . 
Now, this result is far from obvious. It is therefore highly 
significant that when we determine the spectral types of 
the ionising stars in these galaxies we find a value which is 
consistent with this result: the dominant ionising stars are 
early-type B or late-type O which have masses of the order 
of 10 to 20 Mq (Panagia 1973), in good agreement with the 
upper mass cutoff suggested by the model. We have there- 
fore to conclude that the solution suggested by the synthetic 
starburst model is fully consistent with our data. 

But how could star formation be continuous? In order 
to sustain star formation in starburst galaxies, some mech- 
anism other than interaction is necessary. Observations of 
nearby starburst galaxies suggests that most (~ 70%) of 
these galaxies are relatively isolated and cannot be associ- 
ated to on-going interaction events (Telles & Terlevich 1995; 
Coziol et al. 1997b) . In our sample, only five galaxies have a 
peculiar morphology suggesting some sort of direct interac- 
tion. The best example is NGC 1089. The galaxy II Zw 40, 
on the other hand, which also has the youngest burst, can- 
not be easily connected to an interaction event. The same 
is true for most of the galaxies in our sample. More than 
two thirds (19) of our galaxies are isolated early-type spi- 
rals (earlier than Sbc) or galaxies with a compact appear- 
ance (14 galaxies). The rest (10 galaxies) looks like normal 
late-type spirals (Sbc or later). There is, consequently, no 
obvious dynamical cause for the bursts observed in most of 
these galaxies, which suggests that some "internal" regu- 
lating or self-sustained star formation mechanism must be 
active. 

In order for star formation to be continuous over a sig- 
nificantly long period of time, star formation rates cannot 
be very high either. This may be in contradiction with some 
observations in very intense starburst galaxies. This is why 
a sequence of bursts may be a more interesting alternative. 
If the bursts are intense but do not last very long, the reser- 
voir of gas will not be exhaust, allowing ignition of new 
bursts later on. Assuming some self-regulated mechanism, 
this process could probably extend over a relatively long pe- 
riod of time. One can already find in the literature different 
mechanisms capable of producing bursts in sequence, either 
by regulation or propagation of star formation from super- 
novae feedback (Gerola, Seiden & Schulman 1980; Kriigel & 
Tutukov 1993) or multiple small merger events on a few Gyr 
period (Tinsley & Larson 1979) . These two alternatives were 
already found to be consistent with the chemical properties 
of nearby starburst galaxies (see Coziol et al. 1997 and 1999 
for example). 

But what evidence do we have in favour of a sequence 
of bursts? In our spectra we have found some variation of 
burst characteristics with the aperture which suggests the 
bursts are not homogeneous in space. The fact that we do 
not observe systematic trends for these variations suggests 
that their origin may be related to the particular star for- 
mation history of each galaxy. This could be one signature 
of internal processes. Considering that these variations are 
observed over regions a few kpc wide, it is probably more 



realistic to assume also that different star formation events 
occurred during different time intervals (this would be the 
case for star formation propagation). The dispersion of the 
K-band absolute magnitude vs. EW(Br7) relation is highly 
suggestive of some intermittent events, superimposed, per- 
haps, on a more global behaviour. 

However, the strongest evidence in favour of a sequence 
of bursts may be the predominance of early-type B stars 
in the ionised regions of these galaxies. Such stars, with 
masses of the order of 10 to 20 Mq, have lifetime of a few 
10 7 yrs only. For one burst population, it is therefore diffi- 
cult to understand how this particular composition of stars 
is possible. This cannot be explained solely in terms of an 
age effect. Indeed, the main sequence lifetime of these stars 
is only a factor 10 older than those of more massive stars 
and a factor 10 younger than the age of the bursts deduced 
from the synthetic starburst model, which would mean that 
we observe all these galaxies in a very particular and short 
phase of their evolution. Considering the generality of our 
observations, this makes it a very improbable event^|. 

A sequence of short bursts, on the other hand, with- 
out any limit on the upper mass cutoff and extending over 
a relatively long period of time (a few 10 s or 1 Gyr) in 
order to mimic continuous star formation, may reproduce 
such phenomenon. Short bursts would produce few massive 
stars, which would then evolve rapidly and all disappear be- 
fore another burst begin. Depending on the frequency of the 
bursts, the average population would then evolve statisti- 
cally towards the predominance of early-type B stars. Note 
that it is probably impossible to reproduce such stellar pop- 
ulation with constant star formation, unless the IMF upper 
mass cutoff is truncated towards low mass stars. 

To finish, we have concluded that we cannot identify 
what is the mechanism responsible for the H2 emission in 
our galaxies. However, the different time scales over which 
these processes culminate should help us estimate which is 
the more probable one. Supernovae shocks remnants are pre- 
dominant only when a high number of massive stars are pro- 
duced and the bursts are very young (a few 10 4 years). In- 
frared fluorescence, on the other hand, is produced by early 
B stars, which radiate ample photons in the Lyman-Werner 
band (A = 912 — 1108 A). These stars have longer life times 
and should therefore prevail only when the burst is old or, 
for some reason, they are predominant in a burst population. 
The fact that B stars do seem to be over abundant in our 
galaxies suggests that fluorescence must be the main mecha- 
nism responsible for the H2 emission. One would then expect 
to see an increase of fluorescent emission in regions where 
B stars dominate over O stars (Tanaka, Hasegawa & Gatley 
1991). This would explain some of the variation we observed 
when we change the aperture, which would constitute, con- 
sequently, another evidence in favour of the multiple bursts 
scenario hypothesis. 



t It is doubtful that we can explain such rare events based on 
the obvious bias that these are starburst galaxies, considering 
that such galaxies form between 25 to 30% of all the galaxies in 
the nearby Universe. 
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6 CONCLUSION 

Although the solution proposed by the application of the 
synthetic starburst models to our sample of starburst galax- 
ies may look counter intuitive, we have shown that it seems 
to be fully consistent with our data. Our analysis confirms, 
therefore, the results previously obtained by Deutsch & Will- 
ner 1986, Doyon, Puxley & Joseph (1992), Coziol & Demers 
1995, Coziol (1996) and Goldader et al. (1997b). 

However, the solution proposed by the models may not 
represent strong constraints on the duration of the burst or 
the form of the IMF. This is because the models cannot dis- 
tinguish between constant star formation or a sequence of 
short bursts over an extended period of time. This alterna- 
tive, the sequence of short bursts, will have to be tested fully 
before we can conclude on the values of these two important 
parameters in starburst galaxies. 

What ever is the solution, our observations are clear 
on one thing, which is that starburst must be a sustained 
or self-sustained phenomenon: either star formation is con- 
tinuous in time or different bursts happen in sequence over 
a relatively long period of time. The generality of our ob- 
servations suggests that this is a characteristic of starburst 
galaxies in the nearby Universe. 
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